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Abstract 

A layer-by-layer description of the asymmetric lattice gas model for 1//- noise suggested 
by Jensen [Phys. Rev. Lett. 64, 3103 (1990)] is presented. The power spectra of the 
lattice layers in the direction perpendicular to the particle flux is studied in order to 
understand how the white noise at the input boundary evolves, on the average, into 
1 / /-noise for the system. The effects of high boundary drive and uniform driving force 
on the power spectrum of the total number of diffusing particles are considered. In the 
case of nearest-neighbor particle interactions, high statistics simulation results show that 
the power spectra of single lattice layers are characterized by different (3 X exponents such 
that (3 X — > 1.9 as one approaches the outer boundary. 

PACS numbers: 05.40.+j; 74.40+k; 64.60.Ht; 73.50.Td 
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Lattice gas models play an important role in exploring self-organized criticality (SOC) 
[]J which is an attempt to find a general unifying principle explaining the ubiquity of 
1 //-noise and fractality in nature. The idea behind these studies is that l//-spectra 
occur when many-body dissipative systems evolve naturally (without a fine-tuning of 
some parameter) towards a set of states without a characteristic length and/or time 
scale [0, ||, [|, [|, ||. The basic feature of / _/3 -noise, with (3 ~ 1, is the long-range and 
self-similarity of the temporal fluctuations of some physical quantity such as the total 
number of diffusing particles N(t). 

In particular, the asymmetric lattice gas (ALG) model introduced by Jensen 
allows to describe 1 / /-noise in a highly interacting system of particles following diffusive 
dynamics. The ALG model can be considered as a transformer which produces a highly 
correlated signal out of a white noise (or completely uncorrelated signal) at the input 
boundary. It has a simple physical interpretation as a model of transport phenomena in 
type -II superconductors. 

In an attempt to understand how the SOC phenomena arise in lattice gas models, 
such as the ALG model, it is also necessary to carry out a more detailed analysis in 
terms of the power spectra of the individual layers in the model. This is the main goal 
of this report. 

We present here the layer-by-layer description of the Jensen's ALG model for 1//- 
noise by studying the power spectra of the lattice layers in the direction perpendicular 
to the particle flux. We studied how the white noise at the input boundary evolves 
throughout the lattice, characterized by different exponents /3 X , such that the critical 
exponent (3 for the total number of particles (and large times) is close to 1 as first found 
by Jensen || . The effects of high boundary drive and uniform driving force on the power 
spectrum of the total number of diffusing particles are also considered. 

The lattice gas cellular automaton model of Jensen 0] is a system of particles residing 
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on a square lattice A of N x x N y lattice sites. The particles obey an exclusion principle, 
so that there is no more than one particle on a given site at time t. With every site 
r G A a variable n(f; t) taking values 0, 1 is associated. As usual n(r; t) = 1 means that 
there is a particle on site r at time t and n(r; t) = that the site is empty. The particles 
on nearest neighbor sites are repelling each other with a central force of strength J, i. e. 

q 

F int (n(f; t)) = - Jn(r; t) ^ n(f + ef, t)ei , (1) 

i=i 

where e*j , i = 1, . . . , q are the unit vectors to the nearest neighbors (in the calculations 
we put J = 1). In the more general case an additional driving force Fd r {r) could be 
applied to the particles so that the total force acting on a particle is 

Ftatif) = F int {r) + F dT {r) . (2) 

The system evolves at discrete time steps following diffusive deterministic dynamics. The 
configuration {n{f;t)}^ k of the system is updated simultaneously. Every particle I is 
moved one lattice site in the direction of the resulting force, so that the new coordinates 
of the particle are 

nit + i) = n(t) + dt(t) , (3) 

where d is the displacement vector, determined by 



dn 



F x 



dy 



Fy 



(4) 



Here F = \F\ and the square brackets [.] mean taking the nearest integer to the number 
enclosed in the brackets. Equation (3) in combination with (4) implies that diagonal 
moves are also allowed. 

To ensure single occupancy of a lattice site, some additional rules, termed by Jensen 
blocking mechanisms, are applied. These are: 

1) If the site to which a particle attempts to move is occupied, the particle is not moved; 
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2) If there are two particles which attempt to move to the same lattice site, the particle 
to which a larger force is applied wins ; if the forces are equal - neither particle 
moves. 

The boundary conditions are asymmetric. It is assumed that at the left side of 
the system (i.e. at (0, y)) there is a layer of fixed particles . The role of this layer is 
to push the particles in the first column (l,y) into the system. At every time step the 
particles in the first column are first removed and then new particles are introduced with 
a probability p, called boundary drive. In this work, the effect of the lattice boundary 
drive p is studied over a large range of possible values (0-1). The particles can freely leave 
the system over the right edge. Periodic boundary conditions are imposed in y-direction. 
This leads to a net flux of particles through the system. 

The power spectrum Sjv(/) of the total number of particles in the system is defined 
in a standard way as [[7]. §] 

S N (f)=\im |(|X(/,T)| 2 ), (5) 

where X(f, T) is the finite Fourier transform of N(t) (from sample sequencies iV(tj) , i = 
1, . . . ,T) given by 

X(f,T) = f N{t)e- l27rft dt (6) 
Jo 

and the angular brackets <> denote averaging over many realizations of the process. 

If the process is stationary 5jv(.f) can be expressed in terms of the autocorrela- 
tion function of the time signal Rn{t) = lim^oo k J T (N(t)N(t + r)) dt (the so-called 
Wiener-Khintchin's theorem) 

S N (f) = 4 / R N (r) cos (2tt fr) dr . (7) 
Jo 

We focus now on the results obtained for the system from the present simulations. 
Our simulations confirm that the power spectrum has a power low dependence on the 
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frequency in a wide interval of frequencies, i.e. S^{f) ~ 1//^ , with [3 roughly equal to 
1. As can be seen from Fig.l, at small frequencies / < f c , there is a crossover to a white 
noise (Sjsr(f) ~ const.) due to the finite system size. Indeed, f c decreases with increasing 
the system size as ~ 1/N X 2 as predicted in || but with different proportionality constant. 

Having in mind that when the system size N x — > oo the power spectrum is expected 
to be ~ 1/ f 13 , we approximate the power spectrum for the finite system by 

m = / A , (8) 

(l + (f/fc)J 

where A, f c and f3 are fitting parameters. We find that the best fit to our results for 
several system sizes N x x 8, N x — 10, 20, 30, 40, 50, 100 for A, f c and (3 gives: 



A « 0.604/ c _/3 , (3 « 1.09 (9) 
as one would expect from Eq.(8). 

When we change the boundary drive p between and 1, (3 is found to vary slightly. 
However, the limits p — > and p — > 1 need special attention. For the case p — > 0, as 
studied in Refs. [0, 0, it was reported a transition towards f3 ~ 2. Whereas for a high 
boundary drive, e.g. p = 0.99, we have found (3 ~ 1.4 as compared to (3 « 1.05 when 
p = 0.90 for the same system size. Our simulation results also show that for a high 
boundary drive there exists a tendency of (3 to grow as the system size is increased. For 
p = 0.9 and size 60 x 48 we find (3 ~ 1.4 as compared to (3 « 1.05 for the small 10 x 8 
system. 

In the limiting case p = 1 the system is completely deterministic and can be studied 
analytically. The power spectrum is changed to a new type, consisting only of a few 
characteristic frequencies. For a system of length N x = 8, the spectrum consists of 
the following set of frequencies {0, 0.125, 0.250, 0.375, 0.5}. The presence of even a 
very small stochastic element in the boundary drive, i.e. p — > 1, but p ^ 1 produces a 
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qualitatively different picture. 

Let us continue studying now the individual layers in the asymmetric lattice gas 



model under consideration. The averaged occupancy of a layer of the lattice < n{x) > 
is defined as 

/ l T Ny \ 



< n (x) >= (^EE n ( x '2/'*i) ) ' x = l,...,N x . (10) 

X 1 i=ly=l I 

At low boundary drive p = 0.2, we obtain a typical 'saw-like' shape of the density profile 
for small lattice sizes. At this value of p the average density of particles in the statistically 
stationary state of the system remains low, p « 0.28. Since the particles are repelling 
each other, states with alternating high occupancy layers and low occupancy layers are 
encountered more often. As the system size is increased the density distribution gets 
more smooth as shown in Fig. 2, but it can not be well approximated by the linear 
distribution with constant gradient as discussed in 0. An increase in the boundary 
drive smoothens the density distribution. We add that the application of a driving force 
lowers the average occupancy in the system as compared to the case without driving 
force. 

At low boundary drive, the application of a small uniform driving force destroys the 
1/ ' fP type of the power specrtum as seen in Fig.3. For Fd r > F cr , where F cr « 1.75, 
the power spectrum is the same as for a noninteracting system (except for the hard-core 
exclusion) as shown in Fig. 4. At higher values of the boundary drive and small driving 
force the power spectrum displays again the 1/ f 13 behaviour in a small frequency interval. 
Therein (3 takes a value {(3 ~ 1.6) larger than in the case without driving force (/3 ~ 1.05), 
as shown in Fig. 5. This behaviour could possibly be explained in the following way. The 
higher boundary drive (p = 0.9) produces higher average density in the system. A higher 
density leads to stronger correlations and damps the effect of the driving force. This 
point needs further elucidation in view of the analytical results obtained from Langevin 
equation fll0 |. 
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At a moderate boundary drive, p = 0.2, the autocorrelation functions and the indi- 
vidual power spectra of the layers are studied for different lattice sizes, N x x 8, where 
N x = 10, 20, 30, 50, 100 and also 60 x 48. The power specrtum of the number of particles 
in a layer, Si^if), has a form similar to that for the total number of particles in the 
system as seen in Fig. 6. The critical frequency f c at which a crossover to white noise 
occurs depends on the layer selected, i.e. f c = f c (x). The width of the frequency interval 
A/, in which 1/ f 13 holds, depends also on the position x. As x — > N x , the crossover 
frequency changes toward higher frequency values, and the interval Af(x) broadens. It 
is interesting to note that even the first layer (immediately after the input boundary 
layer - whose power spectrum corresponds to white noise) already has (3 ~ 0.8, though 
in a very narrow frequency interval. 

We have calculated the autocorrelation functions Rl x { t ) °f the system layers (shown 
in Fig. 7) by averaging NL x (t)NL x (t + r) over one simulation and then averaging over 
many (up to 1000) independent simulations for different system sizes. Rl x (t) of the 
single layers are found to decay much slower than the autocorrelation function of the 



whole system. As discussed in Ref. [pTTj] , finite system effects on Rl x (t) should be taken 
into account to get the proper behaviour for the infinite system. 

In Fig. 8, (3 r data for several system sizes are shown as a function of the normalized 
layer distance within the lattice r = x/N x . We find that the exponent of the last layer 
is (3{r — 1) « 1.9 - a result very close to the well known result for the Wiener process. 
Furthermore we also see that, within numerical accuracy (less than ~ 5% error), all j3 r 
lie on approximately the same curve. Thus, we find that while the individual power 
spectra of the layers are different and characterized with different /3 X , their collective 
behavior produces a power spectrum for the total number of particles in the system 
with p « 1.1 (Fig.l). 

We also studied the effect of applying a nonuniform driving force on the power 
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spectrum of the total number of particles in the ALG. For moderate boundary drive 
(p = 0.2) we applied additional driving force with a constant gradient AF^r = 0.1. We 
find that the effect on the power spectrum is much less pronounced than in the case of 
a uniform driving force. In fact one can still approximate the spectrum by 1//^, but 
now /3 :=a 1.2 for a system of size 10 x 8 as compared to (5 ~ 1.05 for the same system 
without shear. 

In conclusion we summarize our numerical findings. In an attempt to understand 
how SOC arises in simple lattice gas models of many particle systems, we have carried 
out a more detailed analysis of the ALG model in terms of the power spectra of the 
individual layers in direction perpendicular to the particle flux. The effects of high 
boundary drive and uniform driving force on the power spectrum of the total number of 
diffusing particles have been also considered. In particular, our results for the limiting 
case p — > 1 supplement previous work . Our findings reveal some interesting features 
for the spatial variation of the (local) exponents f3 x . The layer after the boundary layer 
with the white noise power spectrum, is found to display an exponent /3(x = 2) « 0.8 in 
a very narrow frequency interval. Whereas, we find that the exponent of the last layer 
is 0(r = 1) « 1.9. 
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Fig.l: Log- log plot of the power spectra «Sjv(/) of the total number of particles 
against frequency / at low boundary drive p = 0.2 and no driving force: curve (1) 
is for a system of size 10 x 8, curve (2) is for a system of size 50 x 8 together with 
the continuous fitting function F(f). 

Fig. 2: The average occupancy of a column < n(x) > versus column number for a 
system 100 x 8 - curve (1) and for a system 50 x 8 - curve (2) at p — 0.2. 

Fig. 3: The effect of a small uniform driving force F dr = 0.1 on the power spectrum 
for a system 10 x 8 at small boundary drive p — 0.2. 

Fig. 4: The effect of a small uniform driving force F dr = 0.1 on the power spectrum 
for a system 10 x 8 of a noninteracting system (hard-core interaction only). 

Fig. 5: Log- log plot of the power spectra of a system 10 x 8 at high boundary drive 
p = 0.9: curve (1) for F dr =0.1(/?« 1.6) and curve (2) for F dr = (/3 « 1.05). 

Fig. 6: Log-log plot of the power spectra of two layers of a system of size 30 x 8 are 
shown : curve (1) is the power spectrum of the middle layer - x — 15 (f3 ~ 0.95) 
and curve (2) of the last layer - x — 30 (/3 ~ 1.9). 

Fig. 7: The normalized autocorrelation functions Rl x { t ) corresponding to the 
power spectra in Fig. 6: curve (1) is the autocorrelation function of the middle 
layer x — 15 and curve (2) is the autocorrelation function of the last layer x = 30. 

Fig. 8: The dependence of the critical exponent f3 r on the relative distance of a 
layer number r at moderate boundary drive p = 0.2 for system sizes N x x 8 with 
N x = 20,30,40,50,100. Curve (1) is for a 20 x 8 system and curve (2) is for a 
100 x 8 system. 
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